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Abstract 

A coupled map lattice whose topology changes at each time step is studied. We show that the transversal dynamics of 
,the synchronization manifold can be analyzed by the introduction of effective dynamical quantities. These quantities are 
defined as weighted averages over all possible topologies. We demonstrate that an ensemble of short time observations 
pan be used to predict the long-term behavior of the lattice. Finally, we point out that it is possible to obtain a lattice 
with constant topology in which the dynamical behavior is asymptotically identical to one of the time- varying topology. 
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,1. Introduction 

In the physical world there are sets whose elements in- . 
tcract with each other. Neurons [1|, individuals [2|, Joseph- , 
son junctions [3|, and computers [Jl are examples of el- , 
ements that compose these sets called networks. Some . 
■properties of a network derive from its topology, which , 
can be different in each of the abovementioned situations, , 
"namely: regular, random, small- world, and scale- free [5|. 
Although the structure in each network is different, a com- , 
mon feature merges in these mathematical systems: the , 
capability of its elements to synchronize. Aiming to im- , 
■prove the understanding on systems that present interac- ^ 
.tions between their elements, in this letter we dealt with ^ 
■networks whose topology varies with time [2, Iwllf . Based , 
jon the synchronization capability we present exact results 
in order to predict the long-term behavior of these net- ^ 
.works. 

The main idea of this work is the statement of effec- , 
,tive quantities for the analysis of the transversal dynamics 
■to the synchronization mainifold S, from which, for in- 
stance, the entire Lyapunov spectrum can be calculated. ' 
An important point here is that these quantities allow the 
construction of a lattice with an static coupling which, 
for synchronization purposes, exhibits behavior consistent 
with the networks with time-varying topology. We say 
that such effective coupling represents the effective dynam- 
ics of the system. We hypothesized that each realization 
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of the original time-varying network, for a long enough 
observation period, can be replaced by an effective sys- 
tem. Such a system is obtained by calculating a weighted 
average over the set of possible states. 

The manuscript has six main sections including the in- 
troduction. Firstly, we shall present the class of coupled 
map lattice, with time-varying topology, that is the sub- 
ject of this work. The next section outlines the mathe- 
matical tools for calculating the effective quantities. This 
will be followed by the derivation of the expression for 
such quantities and the construction of the effective sys- 
tem. The next section will present a application of our 
results in a coupled map lattice whose coupling between 
the elements is semirandomly determined: the coupling 
probability between two elements in a given instant de- 
pends algebraically on their distance in the lattice. This 
coupling form was studied in Refs. 12, llSj as a model of 
a network of neurons. Our conclusions are left to the last 
section. 

2. Coupled map lattice 

We examined a network called coupled map lattice (CML) 
|14| . which is defined in the following way: let g : w — > w 
be an axiom- jA|^ d— dimensional map defined in w C Mf^. 
Let F : f] — > 17 be a lattice with N coupled maps in the 
form g, with il = w^. This lattice reads 



y„+i = F(y„,77,) = G„f(y„ 



JO) 



(AT-l) 



(1) 



defines the 



in which y„ = x^i ' • • • x;, 

state of the m-th site, f '■'"■' (yn) = g(x„ ), and G„ is a 
dN X dN matrix that depends on discrete time n, in an 
explicit way. This is the coupling matrix. 



^A map that satisfies the axiom- A is hyperboUc and mixing [15 
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2. 1 . Topology si 

We shall restrict the lattice to the linear couplings, as ai 

we can see in the equation ([T]), and consider the lattice sa 

topology varying stochastically with time. However, the 84 

following conditions must be observed: ss 

(i) the lattice is a periodic one; se 

(ii) the topology at each time interval is invariant under s? 

translations in the lattice. ss 

Since there are constraints about the topology, one can 
call this form of semirandom coupling. The above condi- 
tions imply G„ is a circulant matrix. Due to periodicity '' 
imposed to the topology, there is no reason to define a pref- '^ 
erential direction (say, to the left or to the right) and we '^ 
also require the topology to be a symmetric one. These 
features do not depend on the probability rule that de- 94 
fines each possible coupling matrix. The condition that 
the phase space il is a direct product of A'^ subspaces a; is 
satisfied if, and only if. 



dW-l 



[G„],^ > 0, 



E [GJ, 



1, 



{yq,w). (2) 



Herein, q,iv = 0,1,- ■ ■ , dN—l. Thus, G„ is also a stochas- 
tic matrix. We imposed also that the system ([T]) represents '6 
the coupling (direct or indirect) between all the N maps. 97 
That is, the state of an arbitrary site eventually influences ^s 
all other elements in the lattice. If not, the system ([Ij can '^ 
be splited in two uncoupled subsystems. In this case, ouri™ 
results are applied separately to each subsystem. Mathe-id 
matically, this means that the condition < [G/tJom < 1 
must be satisfied, for at least one G^, for each w, in a such'"^ 
way that the infinite product of matrices Gk typically re- 
sults in a strictly positive matrix |16l. llTj. For instance. 



if each site is always coupled with its first neighbors, this 
condition is always satisfied. Here, we consider that the 
network topology uniquely defines the coupling matrix G„ , 
i.e. G„ = G„(T„), in which the topology matrix T„ de- 
scribes the connectivity between sites. 

If we suppose that the rules that specify the coupling'™ 
probabilities are well-defined and do not change over time, 104 
we have 105 

l-L"Jo(Ar-r) - 1 if 



[T, 



7i\()r 



i> Pt 



^— , [T„]ojv =1 
is a parameter that' 



V^n\qyj — [-'-nJo(iL.-gmodAr)' 

in which r = 0,1,..., TV', with N' 
[T„]oo, £, is random variable, and pr 
determines the connectivity between sites and r. Weii2 
remark that the connectivity distribution of each site isi" 
uniquely defined by the parameter pr. If ^ is (5— corrclatedi" 
and with uniform probability distribution in [0,1], then 
Pr defines the coupling probability. Therefore, the proba- 
bility of occurrence of a specific coupling matrix depends 
only on entries of the topology matrix: 

N' 
TTfc = n {pr [T,],^ + (1 - Pr)il - [Tfc]o,.)) • (3) 

r=0 



If the conditions (i) and (ii) are fulfilled, there are, at most, 

N + l 

Kpf — 2~2~ distinct matrices. For each pr ^ oi Pr = 
1, Kn is reduced by a factor 2. The evolution of the 
network is then obtained by the selection of G„ from the 
set {Gk}i;^i according the probability tt^. 

The translacional invariance of the lattice topology im- 
plies all sites share the same connectivity degree k„ at 
the same instant n. Along the temporal evolution of the 
lattice, K„ follows a Poissonian distribution, because we 
assume that the p^'s are independent of each other. Dif- 
ferent distributions can be obtained through a rescaling of 
the probabilities tt^ following their respective connectivi- 
ties. 



3. Methods 

3.1. Chaos synchronization 

We are interested in studying the CML ([T]) from the 
point of view of the synchronization manifold stability. 
Such manifold is defined by 



>S=-^yeri,seaj 



(m) 



sVm 



}■ 



that is, we are considering the complete amplitude synchro- 
nization - all sites share the same state at same time. We 
say that a trajectory {yt}"^o achieves S at time n when, 
for any (5 > 0, we observe dist(y„,5) < S. The mapping 
F keeps S invariant, and the dynamic in this subspace is 
determined by g. 

3.2. Linear stability analysis 

The stability of the synchronization manifold can be 
determined through the Lyapunov exponents' spectrum of 
S, calculated for a typical trajectory in this manifold. The 
m— th Lyapunov exponent is defined by 



A„(so)= lim iln||DF(")(so)e„||, 



(4) 



in which D 



-^ is the dTV— dimensional Jacobian matrix, 
Bm is a typical unitary vector in E'"\E"'+^, and E° = fi 
[l5|. Since the dynamics in S is ruled by the map g, which 
is, by assumption, mixing, Am(so) do not depends on Sg, 
apart from a set of null Lebesgue measure, and we have 
Am(so) = Am. The calculations yielded by Qj give Aq > 
Ai ^ • • • ^ Ajvd-i- The analysis of the stability of 5 is done 
by the examination of the greatest transversal Lyapunov 
exponent, denoted by Ax, which is calculated in n\S. In 
the case of networks whose topology varies over time, the 
scaling behaviour of the Hajnal diameter of the product of 
coupling matrices enables us to obtain A_l 



18 



Now, let us develop Eq. (|4|) for typical trajectories em- 
bedded in S. For the sake of simplicity, we shall consider 
d = 1. Let J„(yo) be the Jacobian matrix applied in ([1]), 
which is calculated over {yt}t'=o'- 



Jn(yo) = n GtDf(y,). 



t=0 



Q^'^ Ijv, where In is the identity 



In S we have Df(yt) 

N X N matrix. Thus, the calculation of the Lyapunov 

exponents given by Eq. (J4|) yields 



we have 



A„i = lim — In i I e 

71— fOO Jl 



Xirn 






(5) 



115 in which Xu > is the Lyapunov exponent of the map 

116 g{s). 

117 Once Gt is a matrix satisfying the conditions (i), (ii 

118 and Eq. ^, we can show that|j 

n-l 

lim rrGt = H, 



119 with [H]q-uj ~ [H]omi, for all q and w. The values la-^^^ 

120 baled by [H]om depends on the observed coupling matriXj^^ 

121 chronological sequence, though it docs not hold for thc^^^ 

122 H form. H is a stochastic matrix, whose greatest eigen-j^^ 

123 value is r^o ^ 1, which is associated to the eigenvector j^^, 

124 ho = -^ [ 1 • • • 1 ] ■ All the others eigenvalues, as-n? 

125 sociated with the eigenvectors hg (q > 0) are equal to zero,i48 

126 rjq = 0. 

127 4. Results and Discussion 

128 4-^- The effective quantities 

In the limit n -> oo, the Lyapunov vectors of S, which 
are parallel to the cigendirections of the Jacobian matrix, 
are determined by hg. For this reason, the Lyapunov ex- 
ponents are given by 

A™= lim iln||e^^"77„J, 



129 with Gm = h„i in ([S]) (the index of h^ arc chosen to provide 

130 a non- increasing spectrum). 

For 77, — > oo we have e^^" — > oo and r]m — > {m > 0). 
If we write 

Om^iVm)'^\ (6) 



^ Our results remain valid under less restrictive conditions with 
respect to the possible coupling matrices. Theorems in Refs. llallTII 
guarantee that the infinite product of matrices that are indecompos- 
able and aperiodic converges to a matrix whose all rows are identicali^i 
(matrix H). A matrix is indecomposable if it contains no subma-i62 
trix of order s X (A^ — s) composed of zeros, whatever s > 0. If 
there is a finite sequence whose product results in a strictly positive 
matrix, i.e. all matrix elements are positive, the infinite product'^* 
of indecomposable matrices conver ges to the matrix H. This result 
was used by Lu et al. in Refs. [la . [iSll in the study of coupled mapi65 
lattices in which the coupling is time-varying. They demonstrated 
the need for the union of all sets contains a spanning tree, that is, 
any pair of elements must be eventually coupled. We extend the 
applicability of expression 10 for all Lyapunov exponents of the syn- 
chronization subspace, from which we could build an effective system 
(autonomous and deterministic). The restrictions that we imposed 
on the topology of the network allowed us to derive analytically an 
expression for the effective quantities, Eq. ((HJl. 



A„i = lim — In | 

n^oo Tl 



Ma 



= \u+ In I 



(7) 



The quantities A™ and (e^^Om) represent the rn-th ef- 
fective exponent and Lyapunov number, respectively, and 
we call 9rn the m-th effective coupling eigenvalue. For 
n 3> 1 we have nA„i, which describe the time-n average be- 
havior of the system. If rjm approaches zero faster than the 
exponential, the limit n — > oo in Eq. ([6]) results 9^ = 0, 
therefore, Eq. ([7]) gives A,„ = — oo. 

Expression ^ gives a general and strong meaning to 
the effective coupling eigenvalues. The 0m (rn > 0) quan- 
tity is obtained by the decay rate to zero of the eigenvalues 
of the matrix H„ = Ilt=o G*- C)nly one condition was im- 
posed, namely lim„_i.oo H„ = H. In other words, for any 
system in the form ([T|), in which all the possible coupling 
matrices satisfies (i), (ii) and Eq. ^ (see [Appendix A[ ), 
the results © and ([71) are formally valid. Thus, in gen- 
eral, our results suggest that the effective dynamics of a 
network whose topology is time-varying can be obtained 
by the analytical or numerical calculation of 9m- 

Since the conditions (i) and (ii) imply all coupling ma- 
trices are circulant, we can write down the exact expres- 
sion for (|n]). This is possible because all circulant matri- 
ces share the same basis and, consequently, commute with 



each other. For this reason, (see Appendix A I 



ATjN 



n 4? 



(8) 



k=l 



in which the productory extends over all possible coupling 
matrices, and Tm is the m-th eigenvalue of the fc-th cou- 
pling matrix. Here and in what follows, the Tm (and 
the 9m) are not indexed by their magnitudes, but accord- 
ing their respective eigenvectors whose components are 
hq' oc exp (— (27ri)(7w/iV)|20|. Notice that, due symmetry, 
Fm = ^N-m ^'^'^ ^11 ^h*3 ^™ '^ith m > are degenerated. 
If for some k, Tm = 0, then 9m = 0, which implies the sys- 
tem is superstable in the ?Ti-th direction. Strictly speak- 
ing, in this situation, any implementation of the system 
eventually results in superstability in the to— th direction. 
The synchronization manifold is transversely superstable 
if, and only if, 9m = in all transversal directions, i.e. for 
every m > 0. In average, the time required for a typical 
trajectory to experience the effects of superstability is of 
the order of tt^"'^. 

4.2. The effective coupling 

Let diag{6'm} be a diagonal matrix formed by the effec- 
tive coupling eigenvalues (|S]) . Let P be an unitary similar- 
ity transformation formed by the Lyapunov vectors of iS. 
Hence, we can make an effective autonomous deterministic 
system from the definition of an effective coupling matrix: 

G = Pdiag{^,„}P^ 
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Figure 1: (color online) Pictorial representation of G elements, in^^^ 
which the color scale stands for the coupling intensity. Note the 
periodicity and translational invariance of the coupling. This matrix 
has the following parameters: N = 31, e = 0.7 and a = 0.6 (elements^oi 
on the main diagonal arc not shown). 202 



From the point of view of the synchronization capacity we 
show that the systems constructed by {G„} and G in (IT|) 
are equivalent. 

2( 

For the systems under consideration, we have [P]qw = 
exp{—{2'Ki)qw/N)/\/N, therefore 



[G]qn 




N' 

2 \J 9s cos 



2tts(w — q)^ 

N 



(9) 



Once cos(-) is an even function with period 27r and^" 
00 = 1, is straightforward that G is a circulant, stochastic 
and symmetric matrix. For this reason, the effective ma- 
trix satisfies all the suppositions about the network. An 
example is found in Figure [T] 

In general, due to the explicit dependence of the Jaco- 
bian matrix on the points of the trajectory, we may for- 
mally demonstrate the equivalence between both systems 
only in the linear neighborhood of S. When such depen- 
dence does not exist, the Lyapunov vectors of all points in 
phase space are identical and our results are globally valid 
in fl 2l|,|22]. From this, we conjecture that the equivalence™ 



between the two systems, for synchronization purposes, is 
generic. Our results arc also valid for networks in which 
the topology varies periodically in time, replacing tt^ by 
the relative frequency associated with k-th topology. 

230 

4.3. Example 

232 

We shall apply the above results to a system in which, 
at each time instant, two sites are coupled with probability 
Pr ~ r^"; being r the shortest distance between sites and 
r, po = 1, and a > a parameter that quantifies the range 



of the coupling. We assume the same value of the coupling 
strength e for all connections. Therefore, the coupling is 
expressed by 



l^n \ ^njlqu 



l + (l-e" 



}l^n]Sqn,]—, (10) 



in which t^„ = 2 ^j,^j^[T„]or is a normalization factor ac- 
cording to 1^. Since, by construction, po = Pi — Ij 
the total number of different topologies for this system 
IS A AT = 2 2 . 

This example describes a network with time-varying 
topology and interactions whose probability of coupling 
between two sites depends on the r distance between them. 
This spatial dependence is recovered by the effective ma- 
trix G in the form of the coupling intensities, as shown in 
Fig.ffl 

We present an analysis on the synchronization of chaos 
in this system, as well as illustrations of the equivalence 
between the original system ([T]) and effective static system. 
In view of this, wc consider the temporal evolution of the 
Euclidean distance between the points of the trajectory 
and S, denoted by d„ = dist(y„,5). As a first examina- 
tion, wc consider the evolution of typical trajectories to S 
for a specific realization of this network. The results for 
the numerical simulations are presented in Figure [H From 
the expressions (|S]) and ([7]) we determine a set of param- 
eters {N,e,a) for which S is transversely stable, in which 
we chose the tent map g{s) = 1 — 2\s — 1/2| {Xjj = In 2) 
to rule the dynamics of each site. The rationale behind 
this choice is the monotonic convergence of the trajecto- 
ries to the transverselly stable synchronized state, when 
static couplings are considered [23|. Thus, we can ensure 
that the irregular time evolution of d„ in Fig. [2]- for typ- 
ical trajectories outside the synchronization manifold - is 
precisely due to the time-dependent coupling. For exam- 
ple, at a given instant n, it is possible that each site is 
only influenced by a few distant sites, which are in very 
different states, causing a departure from S. On the other 
hand, the same mechanism may explain the decrease in 
the synchronization time when the states are very close 
together. Figure [2] also shows the numerical simulations 
for the static performance of the system, in which wc ob- 
serve that the convergence is monotonic. 

We must emphasize that typical trajectories of the sys- 
tem, whose topology is time-varying, present an average 
behavior that is determined by the equivalent static sys- 
tem, for any set of parameters - including the lattice size 
- inside the synchronization domain. This result is shown 
in Fig. [3] in which we analyze the average time required 
for synchronization of a set of typical trajectories. The 
synchronization time increases with the parameter a and 
diverges at its critical value, ac{s,N). For this parameter 
value S becomes transversely unstable. The ac is deter- 
mined by the equation ^ when je'^^^j.l = 1 [2l[, being 



^m>0 



{|^m|}. This critical value □ is indicated by 



^This network can present superstability at e = (1 — ^/N) = e^ 
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Figure 2: (color online) The convergence to S of typical trajectorics26o 
for the time-varying topologies (dashed lines) and for the effectivOjgi 
system (solid lines), with N = 31, e = 0.7 and a = 0.6. 




Figure 3: (color online) The dependence on the range of coupling274 
probability of the average synchronization times, for the time- varying 
system (circles) and the effective static one (diamonds), and its stan- 
dard deviation, given by vertical and horizontal bars, respectively, for^^^ 
iV = 31. The vertical dashed lines are the critical values ac for which277 
S becomes transversaly unstable. 278 



the dashed hnes in Fig. [31 The transition to synchroniza- 
tion, which is characterized by changes in transversal sta- 
bihty of S, occurs at the same point for both systems. In 
this figure, we also analyze the dispersion of synchroniza- 
tion times. The dispersion is computed by the standard 
deviation of several network realizations. This quantity is284 
represented by vertical and horizontal bars in the case of 
systems with varying and static topology, respectively. ^"^ 

286 
287 

5. Conclusions 

288 

In conclusion, we studied here the synchronization 01 
chaos in coupled map lattices whose topology varies semi- 
stochastically with time. Regarding the periodicity, the 
translational invariance of lattice topology, and that alhn 
sites are eventually coupled, we showed that a weighted 
average over all possible topologies determines the dynam- 
ics transverse to the synchronization manifold. This was 



for any a < 00, because global coupling has probability Wgi = 
(A'"'!)"". Typically, this superstability is not observed because the 
set of parameters, which is determined by ess, has null measure in 
parameter space. For e = Ess ± 5e the contribution to the expres- 
sion ^ of the eigenvalue associated with global coupling is given by 
(5e'^ ■' , i.e. it becomes evident only for very small networks. For 
this reason, even for e ?s Ess, no manifestation of this phenomenon 
is observed. 



quantified by the introduction of effective Lyapunov ex- 
ponents, for which we derive an exact analytical expres- 
sion. We also showed that, starting from such effective 
dynamical quantities, a system that is deterministic and 
autonomous - which we call effective system - can be con- 
structed, so that its topology is static. We demonstrated 
the equivalence between the two systems by analyzing the 
stability of the synchronization manifold, and the average 
time needed for the typical trajectories of the network to 
reach it. Although the average behavior is the same for 
both systems, wc showed that the time dependence of the 
coupling can cause great variability in time synchroniza- 
tion for the case of isolated observations. 

Wc found an efficient method to numerically determine 
the effective quantities. This method is based on the eigen- 
value decay rate of the coupling matrices' products, which 
may be a quick way to calculate the Lyapunov spectrum, 
because it usually involves a small number of products of 
matrices (compared to the number of possible topologies) 

Finally, we expect that any network with time- varying 
topology can be represented by a circulant effective matrix, 
provided that the network is periodic and that just the 
probability of interaction (and not the topology) between 
elements is invariant over translations. Consequently, know- 
ing the Lyapunov exponents that characterize the transver- 
sal dynamics to the synchronized state, we can find a good 
estimate for the effective coupling between the elements of 
the system under study. This is possible with the use of 
equation ^ after inversion of equation (O . Note that the 
Lyapunov exponents can be obtained from exact calcula- 
tion, numerical simulations or the analysis of time series 
from experiments. 
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Appendix A. Derivation of Eq. ([8]) 

The restrictions on the topology and coupling restric- 
tions allow an analytical derivation of the effective quan- 
tities. Since any matrix satisfying the conditions (i), (ii) 
and Eq. ([2|) is circulant [20[. A matrix A^ x iV is circulant 
if its elements satisfy 



[G]q,u, — [G]o. 



{w~q mod A^) 



yq,w = 0,l,--- ,N ~1. 



The TO— th eigenvalue of a circulant matrix that is asso- 
ciated to the eigenvector with components hq" ex exp {—{27Ti)qw/N 



is given by 



Af-l 
111=0 



011)6 



~{27Ti)raw /N 



292 If the matrix is circulant and symmetric, this expres 

293 sion reads 



[G]c 



2 5^[GK 



/ 2Trmw\ 



V N 



■ 



in which N' ={N -l)/2. 

Assume that H„ represents a specific sequence 



H„ 



Q{kn-l}Q{k 



. Q{ki} Q{ko) 



328 
329 
330 
331 
332 
333 
334 
335 
336 
337 
338 
339 

(A.I)3" 
341 
342 
343 
344 
345 
346 
347 
348 
349 
350 
351 



(A.2) 



being fc„ a label that identifies a coupling matrix observed, 
at time n, in the set of the Kjy possible coupling matrices. 352 
Since all the matrices have the same basis, they commute^ss 
among themselves and we can rewrite the Eq. (|A.2p in the 
following way: 



g354 
355 
356 
357 

fc=l '™ 

361 

in which i/^ represents how many times the fc-th matrix362 
G*^'^^ was observed, so ^;, Vk = n. In the limit n — > 00 weses 
have lyk = niTk, therefore ^'"' 

365 



Kf 



ptH„p = p^ n (^ 



iWf'p 



fc=i 



KJV 

n (diaglrit^} 



fc=i 



304 Since PtH„P = diaglr;™}, Eq. ([5]) is obtained by 

305 taking the n-th root of expression above. 
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